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Abstract 

In this paper, we derive and analyze a compartmental model for the control of ar¬ 
boviral diseases which takes into account an imperfect vaccine combined with individual 
protection and some vector control strategies already studied in the literature. After the 
formulation of the model, a qualitative study based on stability analysis and bifurcation 
theory reveals that the phenomenon of backward bifurcation may occur. The stable 
disease-free equilibrium of the model coexists with a stable endemic equilibrium when 
the reproduction number, Rq, is less than unity. Using Lyapunov function theory, we 
prove that the trivial equilibrium is globally asymptotically stable; When the disease- 
induced death is not considered, or/and, when the standard incidence is replaced by the 
mass action incidence, the backward bifurcation does not occur. Under a certain condi¬ 
tion, we establish the global asymptotic stability of the disease-free equilibrium of the 
full model. Through sensitivity analysis, we determine the relative importance of model 
parameters for disease transmission. Numerical simulations show that the combination 
of several control mechanisms would significantly reduce the spread of the disease, if we 
maintain the level of each control high, and this, over a long period. 
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1 Introduction 


Arboviral diseases are affections transmitted by hematophagons arthropods. There 
are cnrrently 534 virnses registered in the International Catalog of Arbovirnses 
and 25% of them have cansed docnmented illness in hnman popnlations [1, 2, 3]. 
Examples of those kinds of diseases are dengne, yellow fever, Saint Lonis fever, 
encephalitis. West Nile fever and chiknngnnya. A wide range of arbovirns diseases 
are transmitted by mosqnito bites and constitnte a pnblic health emergency of 
international concern. According to WHO, dengne, cansed by any of fonr closely- 
related virns serotypes (DEN-1-4) of the genns Flavivirns, canses 50-100 million 
infections worldwide every year, and the majority of patients worldwide are chil¬ 
dren aged 9 to 16 years [4, 5, 6]. The dynamics of arboviral diseases like dengne or 
chiknngnnya are inflnenced by many factors snch as hnman and mosqnito behav¬ 
ior, the virus itself, as well as the environment which directly or indirectly affects 
all the present mechanisms of control. 

For all mentioned diseases, only yellow fever has a licensed vaccine. Neverthe¬ 
less, considerable efforts are made to obtain the vaccines for other diseases. In 
the case of Dengue for example, the scientists of french laboratory SANOFI have 
conducted different tries in Latin America and Asia. Thus, a tetravalent vaccine 
could be quickly set up in the coming months. But in any case, it is clear that 
this vaccine will be imperfect. However, the tries in Latin America have shown 
that vaccine efficacy was 64.7%. Serotype-specihc vaccine efficacy was 50.3% for 
serotype 1, 42.3% for serotype 2, 74.0% for serotype 3, and 77.7% for serotype 4 [7]. 
The tries in Asia have shown that efficacy was 30.2%, and differed by serotype [8]. 

Host-vector models for arboviral diseases transmission were proposed in [9, 10, 
11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25] with the focus on the 
construction of the basic reproductive ratio and related stability analysis of the 
disease free and endemic equilibria. Some of these works in the literature focus on 
modeling the spread of arboviral diseases and its control using some mechanism 
of control like imperfect vaccines [19, 20] and other control tools like individual 
protection and vector control strategy [9, 10, 15, 21, 23, 24]. 

In [15], Dumont and Chiroleu proposed a compartmental model to study the 
impact of vector control methods used to contain or stop the epidemic of Chikun- 
gunya of 2006 in Reunion island. Moulay et al. [23] study an optimal control 
based on protection and vector control strategies to hght against Chiknngnnya. In 
[20], Rodrigues et al. simulate an hypothetical vaccine as an extra protection to 
the human population against epidemics of Dengue, using the optimal control. In 
these models [15, 20, 23], 
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(i) the population is constant, 

(ii) the disease-induced death in humans is not considered, 

(iii) the complete stage progression of development of vectors is not considered, 

(iv) none of the above mentioned models takes into account the combination of 
the mechanisms of control already studied in the literature, such as vacci¬ 
nation, individual protection and vector control strategies (destruction of 
breeding site, eggs and larvae reduction). 

The aim of this work is to propose and study a arboviral disease control model 
which takes into account human immigration, disease-induced mortality in human 
communities, the complete stage structured model for vectors and a combination 
of human vaccination, individual protection and vector control strategies to hght 
against the spread of these kind of diseases. 

We start with the formulation of the model, which is an extension of the previ¬ 
ous model study in [26]. We include the complete stage progression of development 
of vectors, the waning vaccine, and four other continuous controls (individual pro¬ 
tection, using adulticides, the mechanical control. Eggs and larvae reduction). We 
compute the net reproductive number A/", as well as the basic reproduction num¬ 
ber, Rq, and investigate the existence and stability of equilibria. We prove that 
the trivial equilibrium is globally asymptotically stable whenever W < 1. When 
J\f > 1 and Ro < 1, we prove that the system exhibit the backward bifurcation 
phenomenon. The implication of this occurrence is that the classical epidemio¬ 
logical requirement for effective eradication of the disease, < 1, is no longer 
sufficient, even though necessary. However considering two situations: the model 
without vaccination and the model with mass incidence rates, we prove that the 
disease-induced death and the standard incidence functions, respectively, are the 
main causes of the occurrence of backward bifurcation. We found that the disease- 
free equilibrium is globally asymptotically stable under certain condition. Through 
local and global sensitivity analysis, we determine the relative importance parame¬ 
ters of the model on the disease transmission. By using the pulse control technique 
in numerical simulations, we evaluate the impact of different controls combinations 
on the decrease of the spread of these diseases. 

The paper is organized as follows. In Section 2 we present the transmission 
model and in Section 3 we carry out some analysis by determining important 
thresholds such as the net reproductive number J\f and the basic reproduction 
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number R^, and different equilibria of the model. We then demonstrate the sta¬ 
bility of equilibria and carry out bifurcation analysis. In section 4, both local and 
global sensitivity analysis are used to assess the important parameters in the spread 
of the diseases. Section 5 is devoted to numerical simulations and discussion. A 
conclusion rounds up the paper. 


2 The formulation of the model 


The model we propose here is based on the modelling approach given in [26, 15, 
16, 17, 18, 19, 23, 24]. It is assumed that the human and vector populations 
are divided into compartments described by time-dependent state variables. The 
compartments in which the populations are divided are the following ones: 

-For humans, we consider susceptible (denoted by Sh), vaccinated (14), exposed 
(Eh), infectious {R) and resistant or immune (i?h); So that, Nh = Sh + Vh + Eh + 
Ih + Rh- Following Garba et al. [19] and Rodrigues et al. [20], we assume that the 
immunity, obtained by the vaccination process, is temporary. So, the immunity 
has the waning rate u. The recruitment in human population is at the constant 
rate Ah, and newly recruited individuals enter the susceptible compartment Sh- 
Are concerned by recruitment people that are totally naive from the disease, and 
immune people whose immunity is lost. Each individual human compartment 
goes out from the dynamics at natural mortality rates fih- The human susceptible 
population is decreased following infection, which can be acquired via effective 

Oif3hv{j]vE-i, -\- A) 


contact with an exposed or infectious vector at a rate \h = 


Nh 


|191 


where a is the biting rate per susceptible vector, jShv is the transmission probability 
from an infected vector (E^ or R) to a susceptible human (Sh)- The probability 
that a vector chooses a particular human or other source of blood to bite can 

1 N-u 

be assumed as —. Thus, a human receives in average ajE bites per unit of 

times. Then, the infection rate per susceptible human is given ^ . 

Nh Ny 

In expression of Xh, the modihcation parameter 0 < < 1 accounts for the 

assumed reduction in transmissibility of exposed mosquitoes relative to infectious 
mosquitoes [19] (see the references therein for the specihc sources). Latent humans 
(Eh) become infectious (R) at rate ■jh- Infectious humans recover at a constant 
rate, a or dies as consequence of infection, at a disease-induced death rate 6. 
Immune humans retain their immunity for life. 

- Following [23], the stage structured model is used to describe the vector 
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population dynamics, which consists of three main stages: embryonic (E), larvae 
(L) and pupae (P). Even if eggs (E) and immature stages (L and P) are both 
aquatic, it is important to dissociate them because, for optimal control point of 
view, drying the breeding sites does not kill eggs, but only larvae and pupae. 
Moreover, chemical interventions on the breeding sites has impact on the larvae 
population (as such as pupae), but not on the eggs [23]. The number of laid eggs 
is assumed proportional to the number of females. The system of stage structured 
model of aquatic phase development of vector is given by (see [23] for details) 



( 1 ) 


P =lL-{9 + fip)P 


V 


Unlike the authors of [23], we take into account the pupal stage in the development 


of the vector. This is justihed by the fact that they do not feed during this 


transitional stage of development, as they transform from larvae to adults. So, the 
control mechanisms can not be applied to them. 

A rate, 6, of pupae become female Adults. Each individual vector compartment 
goes out from the dynamics at natural mortality rates The vector susceptible 
population is decreased following infection, which can be acquired via effective 



contact with an exposed or infectious human at a rate 


where f3hv is the transmission probability from an infected human [Eh or Ih) to a 


susceptible vector (5^). Latent vectors (E^) become infectious {1^) at rate y.,,. The 


vector population does not have an immune class, since it is assumed that their 
infectious period ends with their death [17]. 

Then, we add new terms in the model to assess the different control tools 
studied: 

(i) tti represents the efforts made to protect human from mosquitoes bites. It 
mainly consists to the use of mosquito nets or wearing appropiate clothing 
[24]. Thus we modify the infection term as follows: 


A)j = (1 — q:i)A/i, and A)] = (1 — ai)A„,with 0 < ai < 1; (2) 


(ii) rji and ^2 are eggs and larvae mortality rates induced by chemical intervention 
respectively. 
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(iii) Cm is the additional mortality rate due to the adulticide, 


(iv) 02 is the parameter associated with the efficacy of the mechanical control. 


The above assumptions lead to the following non-linear system of ordinary differ¬ 
ential equations 


Vh 

Eh 

4 

Rh 

Sy 

Ey 

4 

E 

L 

. P 


— Kh + + 'C + h/i) Sh 

= — [(1 ~ + ^ + h/i] 

= [Sh + {l-e)Vh]-{fih + lh)Eh 

= jhEh - (/i/i + + ct)4 

^4 ^^hRh 

= OP - X^Sy - {^y + Cm)Sy 

^yPv (hf T Xv T ^m}Ey 

XvEy (p-t) T Cm^Iy 


— 0-b 

= sE 



{Sy + Ey + ly) “ ( S -|- E T Vl) E 

— {I Pl + ^2)L 


= iL-{e + pp)p 


(3) 


It is important to note that no intervention measure is performed to kill the pupae 
for two reasons: the hrst reason is the fact that at this stage, no food is absorbed by 
the insect, so it is impossible to make her ingest a toxic substance; the second reason 
is the fact that products soluble in water deposits by contact are not selective 
mosquito nymphs and act on all the wildlife of the cottage. 

The description of state variables and parameters of model (3) are given in 
Tables 1 and 2-3. 


2.1 Well posedness of the model 

We now show that the system (3) is mathematically well dehned and biologically 
feasible. We write 

4 := ^ + P/i; 4 := ca p/^; 4 := p/^ 7/^; 4 := P/i + 5 + cr; 

4 := s-F pp pi; 4:=^ + Pl + P2; h := 9 + ftp; ks := fiy + Cm] (4) 
4 ■= Pd + + Cm', Kp '■= ot2TE', Kp := a;2rL; tt := 1 — e. 

System (3) can be rewritten in the following way 

dX , ^ 

— = K{X)X + E (5) 
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Figure 1: A compartment model for vector-borne disease with waning vaccine and 
mosquito aquatic development phase. 
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MX) = 


( -K - ki 

e 

A/i 
0 
0 


u 

-TlXl - k2 
nXh 
0 
0 


0 

0 

-h 

Ih 

0 


Ai(X) 

0 


0 

0 

0 

—A;4 

a 


0 \ 
0 
0 
0 

l-^h ) 


0 1 

A4(X) ) 
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Table 1: The state variables of model (3). 



Humans 


Vectors 

Sh 

Susceptible 

E: 

Eggs 

Vh 

Vaccines 

L: 

Larvae 

Eh 

Infected in latent stage 

P: 

Pupae 

4 

Infectious 

Sy: 

Susceptible 

Rh 

Resistant (immune) 

Ey 

4 

Infected in latent stage 
Infectious 


and 


dl2(X) = 


where Agg = fif, 


^ ~{K + ^s) 
Xy 
0 

Age 

0 

V 0 



0 0 0 

— kg 0 0 

'-fv -kg 0 

Age Age —Ag^ 
0 0 Aiog 

0 0 0 



0 e \ 

0 0 

0 0 

0 0 

—Aio 0 

I —kj j 




and 


Aiq = — + kg] and F = (A/*, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0)'^. 

Fl 

Note that A(X) is a Metzler matrix, i.e. a matrix snch that off diagonal terms 
are non negative [30, 31], for all X G M+h Thus, using the fact that -F > 0, 
system (5) is positively invariant in which means that any trajectory of the 
system starting from an initial state in the positive orthant remains forever in 
The right-hand side is Lipschitz continuous: there exists a unique maximal 
solution. 

By adding the hrst four equations of model system (3), it follows that 


Nh{t) = Ah — fihXh — 5Ih < Ah — ^ihXh 


So that 


0 < Nhit) < — + (xhiO) - —) 


OO, 0 < Nhit) < 

I^h 


Thus, at t 










Table 2: Description and baseline values/range of parameters of model (3). 


Parameters 

Description 

Baseline value/range 

Sources 

Afc 

Recruitment rate of humans 

2.5 day“^ 

|19] 


Natural mortality rate 


day ^ 

[19] 


(67x365) 


in humans 


e 

Vaccine coverage 

Variable 


u 

Vaccine waning rate 

Variable 


e 

The vaccine efficacy 

0.61 

[27] 

a 

Average number of bites 

1 day“^ 

[9, 19] 

Phv 

Probability of transmission of 

0.1, 0.75 day-1 

[9, 19] 


infection from an infectious human 
to a susceptible vector 




ih 

Progression rate from to Ih 


h 11 day ' 

[15, 28] 

6 

Disease-induced death rate 

10-3 day-1 

|19] 

a 

Recovery rate for humans 

0.1428 day-1 

[9, 19] 

Vh,Vv 

Modifications parameter 


0,1) 

[19] 

fJjy 

Natural mortality rate of vectors 


nl day-1 

[9, 19] 

Iv 

Progression rate from to 


^,^\day 1 

[15, 28] 

f^vh 

Probability of transmission of 

0.1, 0.75 day-1 

[9, 19] 


infection from an infectious vector 
to a susceptible human 




e 

Maturation rate from pupae 

0.08 day i 

[15, 23, 24] 


to adult 




f^b 

Number of eggs at each deposit 

6 day-1 

[15, 23, 24] 

Te 

Carrying capacity for eggs 

103,10® 

[9, 23] 

Tl 

Carrying capacity for larvae 

5 X 10^,5 X 10® 

[9, 23] 

t^E 

Eggs death rate 

0.2 or 0.4 

[24] 

I^L 

Larvae death rate 

0.2 or 0.4 

[24] 

/ip 

Pupae death rate 

0.4 
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Table 3: Description and baseline values/range of parameters of model (3). 


Parameters 

Description 

Baseline value/range 

Sources 

s 

Transfer rate from eggs to larvae 

0.7 day“^ 

[24] 

1 

Transfer rate from larvae to pupae 

0.5 day“^ 

[23, 29] 

Vl,V2 

Eggs and larvae mortality rates 
induced by chemical intervention 

0.001,0.3 

[24] 

O-l 

Human protection rate 

|0,1) 


a2 

Efficacy of the mechanical control 

(0,1| 

|15] 


Adulticide killing rate 

[0,0.8] 

[15] 


By adding the equations in Sy, Ey and Ey of system (3), it follows that 

Ny{t) = OP - fiyNy 

So that 

0 < Nyit) = —+ ('iV,(0) - —) 

l-Ly \ / 

Thus, at t —)■ oo, 0 < Ny{t) < —since P < 

jjjykj kj 

Therefore, all feasible solutions of model system (3) enter the region: 

V=\ (Sh, Vh, Eh, 4, Rh, Sy, Ey, ly, E, L, P) e iV;, < ^; E < Ke-, L < Kl; 

„ IKl OIKl] 

P < _ -■ Af < _h K 

- 4 ’ MsJ 

3 Mathematical analysis 

3.1 The disease—free equilibria and its stability 

In the absence of disease in the both population (human and Adult vector), 
i.e. = A(5 = 0 (or Eh = R = Ey = R = 0), we obtain two equilibria 
without disease: the trivial equilibrium (equilibrium without vector and disease) 
Sq = (S')), 14 °, 0, 0, 0, 0, 0, 0, 0, 0, 0) and the disease-free equilibrium (equilibrium 
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with vector and without disease) £i = (S')), 0, 0, 0, Ad), 0, 0, E, L, P) with 


^0 _ yo _ 

" ^ihik2 + 0' " ^^h{k2 + 0' 


N?, = 


KEK^k^ke 


fib{KES + keKi) 

p _ KEKik^k^ks (Ad — 1) p _ KEKik^k^k-iks {M — 1) 


E = 


fibO {Kes + kfiKL) 

KEKLk^k^kjks {M - 1) 


libOl {Kes + kfsKL) 


s {uhlKie + k^k^ksKE) ' 
where J\f is the net reproductive number [23, 32, 33] given by 

^bOls 


M = 


k^k^kYks 


(6) 


( 7 ) 


3.1.1 Local stability of disease—free equilibria 

The local asymtotic stability result of equilibria Sq and Si is given in the following. 
Theorem 3.1. Define the basic reproductive number [34, 35] 


Rn — 


lp{l - aifiPhvPvhUhkbke {jh + kiUh) (jv + ksUv) ( tt ^ + ^ 2 ) a2TETL{f^ - 1 ) 


kskikskgUbSh (C + ^2) (ke^L + sTe) 


( 8 ) 


Then, 


(i) if Af < I, the trivial equilibrium Sq is locally asymptotically stable in V; 

(a) if N > 1, the trivial equilibrium is unstable and the desease-free equilibrium 
Si is locally asymptotically stable in V whenever Ro < 1. 


Proof. See Appendix B. □ 

The basic reproduction number of a disease is the average number of secondary 
cases that one infectious individual produces during his infectious period in a 
totally susceptible population. The epidemiological implication of Theorem 3.1 is 
that, in general, when the basic reproduction number, Rq is less than unity, a small 
influx of infectious vectors into the community would not generate large outbreaks, 
and the disease dies out in time (since the DFE is LAS) [19, 34, 35, 36]. However, 
we show in the subsection 3.2 that the disease may still persist even when Rq < 1. 
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3.1.2 Global stabilty of the trivial equilibrium 

The global stability of the trivial equilibrium is given by the following result: 
Theorem 3.2. If J\f < 1, then So is globally asymptotically stable on T>. 


Proof. To prove the global asymptotic stability of the trivial disease-free equilib¬ 
rium So, we use the direct Lyapunov method. To this aim, we set T = X — TE 
with X = {Sh, 14, Eh, Ih, Rh, Sv, E.^, ly, E, L, P)^ and rewrite (3) in the following 
manner 


dY 

dt 


B{Y)Y. 


The global asymptotic stabilty of So is achieved by considering the following Lya¬ 
punov function C{Y) =< g,Y > where g = (1,1,1,1, 1,1,1,1, —, ^ ^ • 

\ /^6 fJjbSl J 

See Appendix C for the details. □ 


3.1.3 Global stabilty of the disease—free equilibrium 

We now turn to the global stabilty of the disease-free equilibrium Si. we prove that 
the disease-free equilibrium Si is globally asymptotically stable under a certain 
threshold condition. To this aim, we use a result obtained by Kamgang and Sallet 
[37], which is an extension of some results given in [35]. Using the property of 
DFE, it is possible to rewrite (3) in the following manner 

( Xs= Ai{X){Xs-XnFE)+Ai2{X)Xi 

\Xi= A2 {X)Xi ^ > 


where Xs is the vector representing the state of different compartments of non 
transmitting individuals {Sh, 14, Rh, Sy, E, L, P) and the vector Xj represents the 
state of compartments of different transmitting individuals {Eh, R, Ey, R). Here, 
we have X 5 = {Sh, W, Rh, Sy, E, L, Pf, Xj = {Eh, R, Ey, Rf, X = {Xs, Xj) and 
Xdfe := = ( 4 °, V^, 0 , 0 , 0 , N^, 0 , 0 , E, L, Pf, 


with 


Ai{X) 


A? A?\ 

A? A^^^ ) ’ 


/ ~iXh + ^i) 
0 

V 0 


a; 0 0 \ 

-{eXI + 12)0 0 

0 —fJLh 0 

0 0 -{Xl + h) j 
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4'^(x) = 


4"^(x) = 


0 0 0 
0 0 0 
0 0 0 
0 0 0 


^12(X) = 


— ^/C5 + Hb 

s(l-^ 


, ^f>(A') = 
s;\ 


Ke 

4 


0 

ab2r]hS^ 

Nh 

0 


0 0 0 /ife — 

0 0 0 0 

0 0 0 0 


a 

ab2S^ 
' Nb 
0 


abir]^S^ 

n 

CtbiTlyTT 

Nh 

0 


yUfe 1 - 


abiS^b 

abiTiVb 


/Ub 1 - 


0 0 0 

0 0 0 
0 0 0 
0 0 0 

0 0 0 

0 0 0 
0 0 0 


—ks 

0 

abiri^{Sb + nVb) 

abi{Sb + FVh) 


Nb 

Nb 

Ih 

—k4^ 

0 

0 

(E)2'f]bS^ 


— kg 

0 

Nb 

Nb 



0 

0 

Iv 

-kg 


A(X) = 


A direct computation shows that the eigenvalues of ^i(X) have negative real parts. 
Thus the system Xs = — ^dfe) is globally asymptotically stable at 

Xdfe- Note also that A 2 {X) is a Metzler matrix. 

We now consider the bounded set Q: 

g = {(A;,, w, Eb, h, Rh, s,, E,, E, L, P) e M'' : Sb < Nb, W < Nb, Eb < Nb, 

h < Nh, Rh < Nh, Nb = Ab/(nb + b) < Nb < Njl = A^/ fib] 

E<Ke]L<Kl]P<^-^]N,<^-I^\ 

Kj K^Ks I 
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Let us recall the following theorem [37] (See [37] for a proof in a more general 
setting). 

Theorem 3.3. Let Q C U = ISJ x The system (3) is of class C^, defined on 
U. If 

(1) Q is positively invariant relative to (9) . 

(2) The system Xs = ~ ^dfe) is Globally asymptotically stable at 

Xdfe- 

(3) For any x E Q, the matrix . 4 , 2 ( 3 ;) is Metzler irreducible. 

(4) There exists a matrix A 2 , which is an upper bound of the set 

M = { 4 . 2 ( 3 ;) G 444(M) : 3 ; G with the property that if A 2 G A4, for any 
X E Q, such that 4 . 2 ( 3 ;) = 42, then x G x {0}. 

(5) The stability modulus of A 2 , a(42) = maxx<=sp{A 2 )'^^W satisfied a(42) < 0. 
Then the DFE is GAS in Q. 


For our model system (3), conditions (1-3) of the theorem 3.3 are satished. An 
upper bound of the set of matrices 44, which is the matrix 42 is given by 


/ 

—kg 

0 

ohipfiSl + 7rl4°) 





Nh 

Nh 


Ih 

—ki 

0 

0 


ab2r]hSy 

ab2Sl 

— kg 

0 


Nh 

Nh 



V 

0 

0 

Iv 

-kg j 


where Nh = 


A/t 

Ah + d) 


To check condition 
To this aim, let 


(5) 


A 

C 


7h -ki J 

ab2VhSy a&2>S'° \ 

Nh Nh 

0 0 / 


in theorem 3.3, we will use the useful lemma [37] in A. 


abiVvAl + ttV^) 


D 


Nh 

0 

-/Cg 0 
Iv -kg 


a6i(^° + vrf4°) 


Nh 

0 
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Clearly, A is a stable Metzler matrix. Then, after some computations, we 
obtain D — CA~^B is a stable Metzler matrix if and only if 


Rc<l 


( 10 ) 


where 


Rr = 


I a?(l - aif PhvPvhk^kQ + k4r]h) (jv + ksrjy) KEKL{k2 + 7rC(A/' - 1) {fih + 


k3k4kskQ^bik2 + OikeKL + KEs)Ah 


^J-h 


( 11 ) 


We claim the following result 

Theorem 3.4. If J\f > 1 and Ro < Rc < I, then the disease-free equilibrium Si is 
globally asymptotically stable in Q. 

Remark 3.1. From (11), we have 


2 {Rh + dY 


Ri = 


rI 


Rq > -^ 0 ) 


showing that R^ is not necessarily an optimal threshold parameter. 


Remark 3.2. Note that in the absence of disease-induced death, i.e. 6 = 0, we 
have Rc = Rq. This suggests that the disease-induced death may be a cause of the 
occurence of the backward bifurcation phenomenon. 

Remark 3.3. The previous results are of utmost importance, because they show 
that if at any time, through appropriate interventions (e.g. destruction of breeding 
sites, massive spraying, individual protection,...), we are able to lower M or Rq 
and Rc to less than 1 for a sufficiently long period, then the disease can disappear 
[ ISj . 

Theorem 3.4 means that for Rq < Rc < 1, the DFE is the unique equilibrium 
(no co-existence with an endemic equilibrium). If i?,, < -Ro < 1; then it is possible 
to have co-existence with endemic equilibria and thus, the occurrence of backward 
bifurcation phenomenon. 

The backward bifurcation phenomenon, in epidemiological systems, indicate 
the possibility of existence of at least one endemic equilibrium when Rq is less 
than unity. Thus, the classical requirement of Rq < 1 is, although necessary, no 
longer sufficient for disease elimination [19, 38, 39, 40]. In some epidemiological 
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models, it has been shown that the backward bifurcation phenomenon is caused by 
factors such as nonlinear incidence (the infection force), disease-induced death or 
imperfect vaccine [19, 40, 41, 42, 43, 44], To conhrm whether or not the backward 
bifurcation phenomenon occurs in this case, one could use the approach developed 
in [35, 43, 45], which is based on the general centre manifold theorem [46]. We will 
explore this method in the next section. 

3.2 Endemic equilibria and bifurcation analysis 

3.2.1 Existence of endemic equilibria 

We turn now to the existence of endemic equilibria. Let us introduce the following 
quantity i?i = Rl\s=o. We proove the following result 

Theorem 3.5. We assume that Af > 1, then 

(i) In the absence of disease-induced death in human population (5 = 0), model 
system (3) have 

1. an unique endemic equilibrium whenever Ri > 1. 

2. no endemic equilibrium otherwise. 

(a) In presence of disease-induced death in human population (5 > 0), model 
system (3) could have 

3. at least one endemic equilibrium whenever Rq > 1. 

f. zero, one or more than one endemic equilibrium whenever Rq < 1. 

Proof. See appendix D. □ 

Note that case 4 of Theorem 3.5 indicate the possibility of existence of at 
least one endemic equilibrium for Rq < 1 and hence the potential occurrence of a 
backward bifurcation phenomenon. 

3.2.2 Backward bifurcation analysis 

In the following, we use the center manifold theory [21, 35, 43, 45] to explore the 
possibility of backward bifurcation in (3). To do so, a bifurcation parameter 
is chosen, by solving for /3hv from Rq = 1, giving 

__ kshkshfXbd^h (C + h) {hKi + sKe) _ 

a?{l-aiYlIyh^hkbkQf^h + kAPh) (7„ + fcs( tt^ + ^ 2 ) - 1)' 

( 12 ) 


16 



Let Jp* denotes the Jacobian of the system (3) evaluated at the DFE {£i ) and 
with (3hv = I3hv Thus, 

(13) 


- 


Ji 

Js 


J 2 

J 4 


where 


Ji — 


/ —ki 

i 
0 
0 

V 0 


OJ 

-k2 

0 

0 

0 


0 

0 

-ks 

Ih 

0 


0 

0 

0 

—ki 

a 


0 

0 

0 

0 


Ji = 


-fJ-h ) 


( -ks 
0 
0 

Ki 
0 
0 


V 


Jo = 


0 - 


«(1 - oi)l^hvVvSl a(l - a)Pl^Sl 




0 - 


a(l - <^)Phv'^VvV^ 


n 

0(1 - a)/ 3 *^TrV° 


V 0 


J3 = 


0 0 - 


0 0 


a(l - 

a(l - o)PhvK° 

K 

K 

0 

0 

0 

0 

a(l - a)l 3 yhr]hS° 

a(l - a)l 3 yhS° 

K 

K 

a(l - a)/ 3 yhVhSS 

a(l - a)KhS° 


0 

-kg 
Iv 

Ki 
0 
0 

0 0 0 
0 0 0 
0 0 0 


0 

0 

-ks. 

Ki 

0 

0 

\ 


0 

0 

0 

-Kg 

Kg 

0 


0 

0 

0 

0 

-Ki 

I 


^ \ 

0 

0 

0 

0 

-fcr j 


K 

0 

0 

0 

0 


K 

0 

0 

0 

0 


0 


with. = S^+ttVI^, Ki = /ih 
sK 


1 - 


E* 


K2 = h + -^Sl Kg = 

it K 


1 - 


L* 


and Ki= fee + 


Kr 


Note that the system (3), with fJhv = has a hyperbolic equilibrium point 
(i.e., the linearized system (3) has a simple eigenvalue with zero real part and 
all other eigenvalues have negative real part). Hence, the center manifold theory 
[46, 47] can be used to analyze the dynamics of the model (3) near jShv = f^hv The 
technique in Castillo-Chavez and Song (2004) [45] entails Ending the left and right 
eigenvectors of the linearized system above as follows. 

The left eigenvector composants of Jp^ , which correspond to the uninfected 
states are zero (see Lemma 3 in [35]). Thus a nonzero composants correspond to 
the infected states. It follows that the matrix Jp* has a left eigenvector given by 
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V = {yi,V 2 , ... where 


ksK 


Vi=V2=V5=Vq = Vs)= Vio = = 0; ^3 = - 

a(l - ai)(3l^H^ 
iVvkg + 7„) 


Vs; 


Vi = 


a(l - ai)/3yhS^{r]^k8 + 'iv) 


kihN^ 


-Vs, V7 = 


kg 


Vs, Vs = Vs> 0. 


The system (3) has a right eigenvector given by w = (tci, W 2 , ■ ■ ■, wn)^ , where 


wii >0, ws> 0, 

fcy Ki9 ks e kg 

I k^ks 'jv ks 'jy 

'yh(JkskgNjl fih ki 

W5 = -TZ -^-^^ 8 ’ "^4 = —W5, Ws = —W 4 , 

a(l - aijfjyhlJ.h'yvb^iirihki + 'Jh) cr 7^ 


W2 = 


a(l - ai)/3l^{r]yks + 'jy 

7^iV°(fci/c2 - ^ 1 :^) 
u a{l - ai)/3*^^S^ 




Wi = —W2 - 
ki 


kiK 


{r]yW7 + ws). 


Theorem 4.1 in Castillo-Chavez and Song [45] is then applied to establish the 
existence of backward bifurcation in (3). To apply such a theorem, it is convenient 
to let fk represent the right-hand side of the k^^ equation of the system (3) and 
let Xk be the state variables whose derivative is given by the equation for 
fc = 1,..., 11. The local bifurcation analysis near the bifurcation point {Phy = (3^^) 
is then determined by the signs of two associated constants, denoted by Ai and 
A 2 , dehned by 


A = 


n 

vkWiWj 

k,i,j=l 


d^fk{0,0) 

dxidxj 


and 


n 

^2 = VkWi 
k,i=l 


d^fk{0,0) 

dxidcj) 


(14) 


with 0 = (]hv—f3hy It is important to note that in /^(O, 0), the hrst zero corresponds 
to the disease-free equilibrium, £ 1 , for the system (3). Since jShy = (d^y is the 
bifurcation parameter, it follows from 0 = (dhy — (d^y that 0 = 0 when fdhy = id^y 
which is the second component in /fc(0,0). 

It follows then, after some algebraic manipulations, that 


= Ti - Ta 
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with 


a(l - ai)/3^^(2i4°wi + 7rS')Jw2) . . 

Ti = ^- - — -{.VvW^ + '«;8)^^3 


W? 

a(l - ai)(3yhSl 


K 


iVhWs + W4)^+ (^r]hW3 + ^W4 


WgV7, 


To — 2 


and 


a(l — ai)/3vhS^ (\ / X 

-- I 1 ^^hW3 + W4)V7 

«(SS + ^VS) 


. i=3 


— 


NS 


(ri„wy + Wa)v3 


Hence, the coefficient > 0 if and only if 


Fi > r2 (15) 

Note that the coefficient A 2 is antomatically positive. Thus, using Theorem 4.1 in 
[45], the following result is established. 

Theorem 3.6. The model (3) exhibits a backward bifurcation at Rq = 1 whenever 
the inequality (15) holds. If the reversed inequality holds, then the bifurcation at 
Rq = 1 is forward. 

The associated bifurcation diagrams are depicted in Figures 2 and 3. Param¬ 
eter values used in figure 2 correspond to those in Table 4, except A/i = 10, e = 1, 
(Ivh = 0.8, r]h = I, T].,, = 1, a = 0.01428, h = 1, oi = 0.001, 02 = 1, Cm = 0.0001, 
TE = 10®, Ti = 50000. In this case the conditions required by Theorem 3.6, are 
satished: Ai = 0.0114 > 0 and A 2 = 1.1393 > 0. 

Parameter values used in hgure 3 correspond to those in Table 4, except Ah = 
10, I3yh = 0.8, r]h = r]y = 0 = 6 = Cm = ai = 0, 02 = I, Te = 10®, Pl = 50000. We 
also have Ai = —2.4223 < 0 and A 2 = 0.8333 > 0. 

The occurrence of the backward bifurcation can be also seen in Figure 4. Here, 
Rq is less than the transcritical bifurcation threshold {Rq = 0.29 < 1), but the 
solution of the model 3 can approach either the endemic equilibrium point or the 
disease-free equilibrium point, depending on the initial condition. 
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Figure 2: The backward bifurcation curves for model system (16) in the {Ro,El), 
and (-Ro) E*) planes. The parameter (3hv is varied in the range [0, 0.2810] to allow 
Rq to vary in the range [0, 1.5]. Two endemic equilibrium points coexist for values 
of Rq in the range (0.2894, 1) (corresponding to the range (0.0105, 0.1249) of 
/3/i^). The notation EE and DEE stand for endemic equilibrium and disease free 
equilibrium, respectively. Solid line represent stable equilibria and dash line stands 
for unstable equilibria. 
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Figure 3: The forward bifurcation curves for model system (16) in the 

and {Rq, E*) planes. Solid line represents stable equilibria and dash line stands for 

unstable equilibria. 


From theorem 3.5, item (i), it follows that the disease-induced death in human 
(6) may be a cause of the occurence of backward bifurcation phenomenon. In the 
following, we show that the backward bifurcation phenomenon is caused by the 
disease-induced death in human and the standard incidence functions (A)j and A))). 
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Figure 4: Solutions of model (3) of the number of infected humans, Eh, and the 
number of infected vectors, Ey, for parameter values given in the bifurcation di¬ 
agram in Figure 2 with f3hv = 0.0105. So Rq = 0.29 < 1, for two different set of 
initial conditions. The hrst set of initial conditions (corresponding to the solid tra¬ 
jectory) is Sh = 700, 14 = 10, Eh = 220, 4 = 100, Rh = 60, Sy = 3000, Ey = 400, 
4 = 120, E = 10000, L = 5000 and P = 3000. The second set of initial condi¬ 
tions (corresponding to the dotted trajectory) is 4 = 489100, 14 = 10, Eh = 220, 
h = 100, Rh = 60, Sy = 3000, Ey = 400, 4 = 120, E = 10000, L = 5000 and 
P = 3000. The solution for initial condition 1 approaches the locally asymptot¬ 
ically stable endemic equilibrium point, while the solution for initial condition 2 
approaches the locally asymptotically stable DFE. 


3.3 The different causes of the backward bifurcation 

The occurrence of backward bifurcation phenomenon in epidemiological models, 
is caused by three factors: the presence of an imperfect vaccine, the presence 
of the death induced by the disease, and the standard incidence rates. In this 
section, we will consider two variants of the model (3) (the corresponding model 
without vaccination, and the corresponding model with mass action incidence), to 
determine the causes of this phenomenon. 
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3.3.1 Analysis of the model without vaccination 

The model without vaccination is given by 


Sh — + hh) Sh 

Eh = — (h/i + lh)Eh 

ih 
Rh 
5. 

E„ 


= 'JhEh - (/i/x + 5 + cr)4 
^4 ^^hRh 

= 9P - \lSy - {jJy + Cm)Sy 
(/^ii T 'y^ T Cyn'jEy 
'IvEy (/^tj T Cm)4 
E \ 

{Sy + Ey + ly) “ (S T E T )-^ 


(16) 


E 


L \ 


E = lib 
L =sE . ^ 

Oi2^L 

P =iL-{e + iip)p 


J 


{I + lip + Tl2)L 


where and A^ are given at section 2. Model system (16) is defined in the 
positively-invariant set 


El — {(‘S'/l, Eh, Ih, Rh, Sy, Ey, ly, E,L,P) G 


Nh < Rh/ h/i) E < Ke', L < Kp] P < 


IKi 

kv 


no 


;Ny < 


OIKl 

kjks 


Without lost of generality, we assume that Af > 1. The corresponding disease-free 
equilibria of model (16) are given by = (A°, 0, 0, 0, 0, 0, 0, 0, 0, 0) and 
gnv _ g^ g^ g^ g^ ^ p'^ witli ^ and E, L and P are the 

same, given by (6). The associated next generation matrices, Ei and Vi, are, re¬ 
spectively, given by 

/ 0 0 a(l - aPPhvVv a(l - ai)(Ihv \ 


Ei = 


0 0 

a(l - aPPyhVvNy a(l - a^flyhN^ 


Vi = 


K 

V 0 

/ ^3 0 0 0 \ 

-7fe ^4 0 0 

0 0 fcg 0 

\ 0 0 -^y kg y 


N/i 

0 


0 

0 

0 


0 

0 

0 


and 


/ 
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It follows that the associated reproduction number for the model without vac¬ 
cination, denoted by Rnv = is given by 


Rn.ij 


' 0^(1 - ai)‘^/3hv/3vh{'lh + k4Vh)hv + hVv)N^ 


hkikshNjl 


(17) 


Using Theorem 2 of [35], we establish the following result: 

Theorem 3.7. Assumed that M > 1. For basic arboviral model without vaccina¬ 
tion, given by (16), the corresponding disease-free eguilibrium is LAS if Rnv < 1? 
and unstable if Rnv > 1- 


Existence of endemic equilibria 

Here, the existence of endemic equilibria of the model (16) will be explored. 
Let us set the following coefficients 


_ {2ks{k3k4 - S'yh) + {r]hk4 + 'yh)afih{,l - (yi)/3vh} 

" ~ ksk4ks 

d2 = -kgUbAhisKE + k^KR) {kgki - SyR {{ghki + "fRaghR - 

di = klkjkgkgisKE keKRfibRhLh{Rlv “ Rc), 
do = klkjkgkgRKE koKRfibAhLl {RL “ 1 ) • 

We claim the following: 


ai)/3vh -t {kgki - SyRkR < 0, 


(18) 


Theorem 3.8. The arboviral diseases model without vaccination (16) has: 

(i) a unigue endemic eguilibrium if do > 1 Rnv > 1; 

(a) a unigue endemic eguilibrium if di > 0, and do = 0 or d\ — 4(72^0 = 0; 

(Hi) two endemic eguilibria if do < 0 (i.e. Rnv < 1), di > 0 (i.e Rf„ > RR and 
d\ — 4(72^0 > 0; 

(iv) no endemic eguilibrium if do <t) (i.e. Rnv < Ij and h = 0. 


(v) no endemic eguilibrium otherwise. 

Proof. Solving the equations in the model (16) in terms of XR and Xf*, gives 


Sl = 




Lh X 


C,* ? 

h 


K = 


\c,* Q* 


T* — 


ihxrsi 

kok4 


ID* _ 

— 


^InXrSl ^ 


(19) 
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and 


OP p. _ OPK* _ 7JPK’* 

" (Ar + fcg)’ hiXr + ksY ^ kshiXT + ksY 
^ ^ fibOKEP ^ ^ _ ^^,esKEKLP _ 

{k^kgKE + iJibOP) ’ k^KLik^ksKE + /ibOP) + s^bOKEP' 

Substituting (19) and (20) into the expression of and A* and simplifying, shows 
that the nonzero equilibria of the model without vaccination satisfy the quadratic 
equation 

d2{^h*y + diXY* + do = 0 (21) 

where di, i = 0,1, 2, are given by (18). 

Clearly, ^2 < 0 and do > 0 (resp. do < 0) if Rnv > 1 (resp. Rnv < !)• Thus 
Theorem 3.8 is etablished. □ 

It is clear that cases (ii) and (hi) of theorem 3.8 indicates the possibility of 
backward bifurcation (where the locally-asymptotically stable DFE co-exists with 
a locally-asymptotically stable endemic equilibrium when Rnv < 1) in the model 
without vaccination (16). 

This is illustrated by simulating the model with the following set of parameter 
values (it should be stated that these parameters are chosen for illustrative purpose 
only, and may not necessarily be realistic epidemiologically): = 5, (3hv = 0.03, 

Vh = Vv = ^, d = I, a = 0.01, Cm = 0.1, Yvh = 0.4, ai = 0.7 and 02 = 0.5. All 
other parameters are as in Table 4. With this set of parameters, Rc = 0.0216 < 1, 
Rnv = 0.2725 < 1 (so that Rc < Rnv < !)• It follows: d2 = —0.0263 < 0, 
di = 4.8763 X10-^ and do = -3.5031 x 10“^ so that d?-4d2do = 2.0093x 10"^ > 0. 
The resulting two endemic equilibria = {S^, E^, R^, S*, E*, /*, E, L, P), are: 

= (281, 70, 5,1207, 5739,182,44, 22180,10201, 9977) which is locally stable and 
= (6333, 67,4,1147, 5936, 37, 2, 22180,10201, 9977) which is unstable. 

The associated bifurcation diagram is depicted in figure 5. This clearly shows 
the co-existence of two locally-asymptotically stable equilibria when Rnv < 1, 
conhrming that the model without vaccination (16) undergoes the phenomenon of 
backward bifurcation too. 

Thus, the following result is established. 

Lemma 3.1. The model without vaccination (16) undergoes backward bifurcation 
when Case (in) of Theorem 3.8 holds. 
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Stable EE 



Figure 5: The backward bifurcation curves for model system (16) in the E^), 
and {Rnv,E*) planes. The parameter is varied in the range [0, 0.9090] to 
allow Rq to vary in the range [0, 1.5]. Two endemic equilibrium points coexist for 
values of Rq in the range (0.2286, 1) (corresponding to the range (0.0211, 0.4040) 
of Phv)- The notation EE and DEE stand for endemic equilibrium and disease free 
equilibrium, respectively. Solid line represent stable equilibria and dash line stands 
for unstable equilibria. 
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Non-existence of endemic eqnilibria for Rnv < 1 and <5 = 0 

Lemma 3.2. The model (16) without disease-induced death (d = t)) has no en¬ 
demic equilibrium when Rnv, 5 =o < 1? o,nd has a unique endemic equilibrium other¬ 
wise. 

Proof. Considering the model (16) withont disease-indnced death in hnman, and 
applying the same procednre, we obtain that the nonzero eqnilibria of the model 
withont vaccination satisfy the linear eqnation 

PiKi* + Po = 0, 

where pi = kQkiQKi 2 apbJdhPh{^ - + k^{ph + a)ksk^Ki 2 PbJdh and 

Po = -pihk 2 .kAksk^Ki 2 Pb^h {RIv,5=o “ l) • 

Clearly, pi > 0 and po > 0 whenever Rnys=o < 1, so that Xf* = —— < 0. 

Pi 

Therefore, the model (16) withont disease-indnced death in hnman, has no endemic 
eqnilibrinm whenever <1. □ 

The above result suggests the impossibility of backward bifurcation in the model 
(16) without disease-induced death, since no endemic equilibrium exists when 
Rnv, 5=0 < 1 (and backward bifurcation requires the presence of at least two endemic 
equilibria when Rnv^s=o < 1) [19, 40]. To completely rule out backward bifurcation 
in model (16), we use the direct Lyapunov method to proove the global stability 
of the DFE. 

Global stability of the DFE of (16) for h = 0 

Dehne the positively-invariant and attracting region 

V2 = {(,S;„ Eb, h, Rh, s,, 4, E, L, P)eV,:Sh< iV^; S, < N^} 

We claim the following result. 

Theorem 3.9. The DFE, Sf", of model (16) without disease-induced death (d = 
t)), is globally asymptotically stable (GAS) in V 2 if Rnv, 5=0 < 1- 

Proof. Consider the Lyapunov function 

Q = qiEb + q2lh + Qs^v + QaIv 
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where 


Qi = -p, (12 = 

h 


0^(1 - aif I5hvl^vh{flvh + lv)N^ _ a(l - ai)l^hv{rivh + lv)N'^ 

-) Qz — - 


Qa = 


(X (1 1) l^hv 


ksk^kgkgN^ 


kskgkgN^ 


ksks 

The derivative of Q is given by 
G = qiEh + Q2ih + QsEv + QaIv 

— Qli^h^h k^Efi^ -|- Q2{'yhEh k4.Ih) T kgE^^ QA{'yvEy k^I-u) 

= g'i((l - ai)\hSh - ksEh) + q2{lhEh - kdh) + g'3((l - ai)\vSv - kgE^) + qAijvEv - ksE) 


l^h 

“’■v: 

, h/i 


tt(l Ot\)l3}iy{TluEy + Iy)S}i k^Eh 

k^h 

ft(l (y.\^l5-ijii{TiiiE}i -j- IfijSu kg Ey 

l^h 


+ q2{7hEh - k4lh) 
E Q4{'yvEy k^Iy) 


= “ o:i)/3hvr]vShEv + k^ail - ai)PhvShIv 

ks Eh ks Ah 

1 , , a‘^{'i--0!if/3hvl3vhirivk8 + 7v)N^ ^ a^{l - ai)^(3hv/3vh{r]vk8 +'yv)E^ 

- hh^ - 

^ ail-aMv.k8 + ,.)NS ^^^^ _ 


kskskgN^ 


Ah 


a{l-ai)f3hv{r]vk8+'yv)N^ Hh \q o t a(\ - ax)^hv{r\vk8 E 

H- , , , ..n -;—a(l - ai)d„/i6,;4- , , , - kgE„ 


kskskgN^ 


Eh 


kMgN^ 


a{l-ai)Phv a{l-ai)/3hv, . 

“T , , '^yh/y j ^ kgly 


fcs/cs 


fcs/cs 


< 


+ 


1.-1 , a{l-ai)l3hv{Vvk8 + 7 v)El° ' 

%“(i - - 


Ev + 


a(l - ai)l3hv - 


a(l - ai)f3hv{r]vk8 + 7v)E^ fih 


k3k8kgN^ 


A; 


a(l - ai)l3vh'nhN^ 


0^(1 - ai)‘^Phv/3vh{'nvk8 + 7v)N^ ^ 1 

k^k^kskgN^ ks " 


Eh 


+ 


a{l - ai)f3hv{r]vk8+7v)N^ fih ^ ^rO 0^(1 - ai)^/3/ii,/3i,/i(??i,A:8 + 7^)^° 


k3k8kgN^ 


E, 


-a(l - ai)/3„/iiV^ - 


k3k4k8kgN? 


k/x 


Ih 


— ik^nv,5=0 ~ 


We have ^ < 0 if Rnv,s=o < 1, with ^ = 0 if 7^i = 1 or i?/i = 0. Whenever Eh = 0, 
we also have = 0, = 0 and E = 0. Substituting Eh = E = = E = ^ 



the first, fourth and hfth equation of Eq. (16) with 5 = 0 gives Sh{t) —?• S'® = iV°, 
Rh{t) 0, and Sv{t) —)■ 5° = iV° as t —)■ cxd. Thus 

[Shit), Eh{t), hit), Rhit), Sht),Eht),hit),Eit), Lit), Pit)] ^ iNl 0,0,0, iVO, 0,0, E, L, P) 
as t —7> oo. 

It follows from the LaSalle’s invariance principle [48, 49, 50] that every solution of 
(16) (when < 1 ), with initial conditions in V 2 converges to Tf'", as t —)■ cxd. 

Hence, the DEE, T™, of model (16) without disease-induced death, is GAS in V 2 

if 7^L,,=o < 1 - □ 

3.3.2 Analysis of the model with mass action incidence 

Consider the model (3) with mass action incidence. Thus, the associated forces of 
infection, \h and h, respectively, reduce to 

^mh ChiVvRv T A) and Xmv CviVhEh T L/j), (22) 

where, Ch = a(l — Oii)l3hv and = a(l — ai)hh- The resulting model (mass 
action model), obtained by using (22) in (3), has the same disease-free equilibria 
given by ( 6 ). Without lost of generality, we consider that Af > 1. The associated 
next generation matrices, and VA are given by 

( h 0 0 0 \ 

- 7/1 A;4 0 0 

0 0 fcg 0 ’ 

\ 0 0 - 7 ^ kg 

where = S^ + ttV^. It follows that the associated reproduction number for the 
mass action model, denoted by Ro,m = piFmVrn^), is given by 

flo,™ = (23) 

where 

_ f “(1 “ 0!i)/3hvRh ilh + kiVh) (ttI + ^ 2 ) ^ AT 7 m _ f “ Cil)Pvh ilv + ksVv) OP 
Phk3k,h + k2) 

Using Theorem 2 of [35], the following result is established: 

Theorem 3.10. Assume that Af > 1. For the arboviral disease model with mass 
action incidence, given by (3) with ( 22 ), the DFE, Si, is LAS if Ro,m < 1; and 
unstable if Ro,m > 1 


/ 


F = 

-*■ m 


0 

0 

0 


0 

0 

0 


ChVvH^ 
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0 

0 


ChH<^ \ 
0 
0 
0 


U™, = 
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Existence of endemic equilibria Solving the equations in the model (16) in 
terms of and gives 


c* _ 

^mh — 


+ ^ 2 ) 


T* = 

mh 


and 


^mh(^2 + vr(fci + A^;)) + kik2 - o;^’ 
P* ^lhX^^hS*mh 


XhXf^ i-’mh 




3M 


^mh 


h'h^3^4 


y:.h = 




mh 




77 '* _ 

-^mh 


\C,* 


(24) 


5* _ 


9P 


(A5)^; + fc8)’ 


E* = 


/3 P \c,* 
'^mv 


kg{Xmv + A;8) ’ 


J* = 




c,* 

mv 


kskg{Xmv + fcs) 


(25) 


Substituting (24) and (25) into the expression of X’^f^ and X^^ and simplifying, 
shows that the nonzero equilibria of the model without vaccination satisfy the 
quadratic equation 

^2(X'^l)^ + eiX'^l + Co = 0 , (26) 

where e^, i = 0 , 1 , 2 , are given by 


62 = kgkgTT [{'fh + ki7]h) CyAh + k^k^ks] 
ksk^kjkgKTr 


ei = 


{Rem Ro,m) y 


(vr^ + k2 

60 = hkiklkgK (1 - i?o „,) , 


with K = kik 2 — ^uj > 0 and 


Rem. 


[(7/4 + hVh) (tt^ + k2)AhCy + (/ciTT + k2)k3k4,k8\ (tt^ + /C 2 ) 

ksk^kgKTT 


62 is always positive and 60 is positive (resp. negative) whenever i?om is less 
(resp. greather) than unity. Thus, the mass action model admits only one endemic 
equilibria whenever i?om > 1 - 

Now, we consider the case i?om < 1- The occurence of backward bifurcation 
phenomenon depend of the sign of coefficient 61 . The coefficent 61 is always pos¬ 
itive if and only if < Rem- It follows that the disease-free equilibrium is 
the unique equilibrium when Af > 1 and Rem < 1- Now if Rem < Rgm < I) 


then in addition to the DFE £ 1 , there exists two endemic equilibria whenever 
Am = ef - 46260 > 0. However, Rem < i?o,m < 1 ^ Rem < I ^ /3vh < 
k3k4ks{^UTT + fc2(vr^ + k2)) 


a(l - CKi) {-fh + kiT^h) (tt^ + fc2)« + k2)Ah 


< 0. Since all parameter of model 
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are nonnegative, we conclude that the condition Rem < -^om < ^ hold. 

And thus, the model with mass-action incidence does not admit endemic equilibria 
for Rem ^ 1* 


Global stability of the DFE for the model with mass action incidence 

Since the DFE of the model with mass action incidence is the unique equilbrium 
whenever the corresponding basic reproduction number R^ ^ is less than unity, it 
remains to show that the DFE is gas. To this aim, we use the direct Lyapunov 
mehod. 

Let us dehne the following positive constants: 



CHH^{r],ks + 'y.)C,S^, 


hh hh ’ 


P3 = PiChH^ 


iVvh + 7^) 

hh 


, Pa = 


ChH^ 
hh ' 


Consider the Lyapunov function 


C =piEh + P2h + PsEy + pjy 


The derivative of C is given by 

t = piEh + P24 + PsK + pJv 

= {piChPvH -h p^'yy - P‘ih)Ey {piChH - Pih)Iy 
+ {psCvPhSy + P2lh - Pih)Eh -f {psCySy - p2kA)h 

Replacing p*, i = 1,.. .4 by their respective term, and using the fact that H = 
(4 + <H^ = {Si + TiV^) and < Nl in 

^3 = {{Sh: Eh, Ih, Rhi Sy, Ey, ly, E, L, P) G V : 

Nh <—,Sy< Nl = 9P,E< Ke,L< Kl,P< , 

Ph kjks J 


we obtain E < - 1) E,. 

We have £ < 0 if Ro,m < 1, with £ = 0 if 4o,m = 1 or = 0. Whenever Eh = 
0, we also have R = 0, Ey = 0 and R = 0. Substituting Eh = R = Ey = R = 0 
in the hrst, fourth and hfth equation of Eq. (16) with mass action incidence gives 
Sh{t) —!• S')), Vh{t) —)• Vh, Rh{t) —!• 0 , and Sy{t) —)■ S')) = Nl as t —)• oo. Thus 

[Sh it) ,Vh{t),Eh{t),R{t),Rh it) , Sy it) , Ey it) ,Iy{t),E{t),L{t),P{t)] 

{Sh^Vh ,0,0,0, NI ,0,0, E, L, P) as t ^ oo. 
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It follows from the LaSalle’s invariance principle [48, 49, 50], that every solution 
of (3) with mass action incidence, with initial conditions in converges to the 
DFE, as f —)■ oo. Hence, the DFE, of the model with mass action incidence, is 
GAS in P 3 if 7^o,m < 1- 

Thus, we claim the following result. 

Theorem 3.11. The DFE, 81 , of the model (3) with mass action incidence, is 
globally asymptotically stable (GAS) in if Ro,m < 1- 

Thus, the substitution of standard incidence with mass action incidence in the 
arboviral model (3) removes the backward bifurcation phenomenon of the model. 
It should be mentioned that a similar situation was reported by Garba et al. in [19] 
and by Sharomi et al. in [40]. 

We summarize the previous analysis of Subsection 3.3 as follows: 

Lemma 3.3. The main causes of occurence of backward bifurcation phenomenon 
in models (3) and (16) are the disease-induced death and the non-linear incidence 
rates. 


4 Sensitivity analysis 

As shown in the previous sections, model (3) may admit single or multiple steady 
states according to the value of the basic reproduction number Rq. In turn, Rq 
depends on the parameters of the model. Usually there are uncertainties in data 
collection and estimated values, as for our model, and therefore it is important to 
assess the robustness of model predictions to parameter values and, in particular, 
to estimate the effect on Rq of varying single parameters. To this aim, we use 
sensitivity analysis and calculate the sensitivity indices of Rq to the parameters in 
the model using both local and global methods. 

4.1 Local sensitivity analysis 

The local sensitivity analysis, based on the normalised sensitivity index of Rq (see 
[51]), is given by 

. _ ORq 
^ Rq aT 

where T denotes the generic parameter of (3). 

This index indicates how sensitive Rq is to changes of parameter T. Glearly, a 
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Table 4: Parameter values using to compute the sensitivity indices of Rq. 


Parameter value Parameter value Parameter value 



0.01 

s 

0.7 

Pvh 

0.75 

t^b 

6 

V2 

0.3 

Te 

10000 

Ap 

0.4 

Ae 

0.2 

Tl 

5000 

e 

0.08 

e 

0.61 

0!2 

0.5 

i 

0.5 

Ah 

2.5 


1 

67*365 

a 

1 

Phv 

0.75 

riv 

0.35 

f^V 

1 

30 

Al 

0.4 

a 

0.1428 

Ih 

1 

14 

Vh 

0.35 

Iv 

1 

21 


0.5 

u 

0.05 

Vi 

0.001 

5 

0.001 

CXl 

0.2 




positive (resp. negative) index indicates that an increase in the parameter value 
results in an increase (resp. decrease) in the Rq value [51]. 

For instance, the computation of the sensitivity index of Rq with respect to a 
is given by 


a dRo 
Rq da 


1 > 0 . 


This shows that Rq is an increasing function of a and the parameter a has an 
influence on the spread of disease. 

We tabulate the indices of the remaining parameters in Table 2 using parameter 
values on Table 4. The results, displayed in Table 5 and Figure 7a. The parameters 
are arranged from most sensitive to least. The model system (3) is most sensitive 
to a, the average number of mosquitoes bites, followed by /r„, e, s, A^, f3hv, f^vh, 
Te, Tl and a2- It is important to note that increasing (decreasing) a by 10% 
increases (decreases) Rq by 10%. However, increasing (decreasing) the parameters 
/i^ by 10% decreases (increases) Rq by 9.190%. The same reasonning can be done 
for other parameters. 


4.2 Uncertainty and global sensitivity analysis 

Local sensitivity analysis assesses the effects of individual parameters at particular 
points in parameter space without taking into account of the combined variability 
resulting from considering all input parameters simultaneously. Here, we perform 
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Table 5: Sensitivity indices of Rq to parameters of model (3), evaluated at the 
baseline parameter values given in Table 4. 


Parameter 

Index 

Parameter 

Index 

Parameter 

Index 

a 

+ 1 

a 

-0.2911 

e 

-0.0566 

fly 

-0.9190 


-0.2757 

UJ 

+0.0565 

e 

-0.6223 


-0.25 

I^E 

-0.0171 

s 

+0.5172 

9h 

+0.2067 

6 

-0.0020 

Ah 

-0.5 

Ih 

-0.2064 

hi 

-0.0000858 

Phv, Pvh, Pe, Pl, 0^2 

+0.5 

Vv 

+0.1207 



f^h 

+0.4996 

Iv 

+0.1174 



/ip 

-0.4810 


-0.1026 



9 

+0.4810 

f^b 

+0.0772 



1 

+0.4489 

V2 

-0.0770 




a global sensitivity analysis to examine the models response to parameter variation 
within a wider range in the parameter space. 

Following the approach by Marino et al. and Wu et al. [52, 53], partial rank 
correlation coefficients (PRCC) between the basic reproduction number Rq and 
each parameter are derived from 5,000 runs of the Latin hypercube sampling (LHS) 
method [54], The parameters are assumed to be random variables with uniform 
distributions with its mean value listed in Table 4. 

With these 5,000 runs of LHS, the derived distribution of Rq is given in Figure 6. 
This sampling shows that the mean of Rq is 2.0642 and the standard deviation is 
2.6865. The probability that i?o > 1 is 54.86%. This implies that for the mean of 
parameter values given in Table 4, we may be conhdent that the model predicts a 
endemic state. 

We now use sensitivity analysis to analyze the influence of each parameter on 
the basic reproductive number. From the previously sampled parameter values, we 
compute the PRCC between Rq and each parameter of model (3). The parameters 
with large PRCC values (> 0.5 or < —0.5) statistically have the most influence 
[53]. The results, displayed in Table 6 and Figure 7 (b), show that the parameters 
«!, the human protection rate, has the highest influence on Rq. This suggests that 
individual protection may potentially be the most effective strategy to reduce Rq. 
The other parameter with an important effect are a 2 , (^hv, f^vh and 9. 

We note that the order of the most important parameters for Rq from the local 
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Figure 6: Sampling distribution of Rq from 5,000 runs of Latin hypercube sampling. 
The mean of Rq is 2.0642 and the standard deviation is 2.6865. Furthermore, 
P(i?o > 1) = 0.5486. 


sensitivity analysis not match those from the global sensitivity analysis, showing 
that the local results are not robust. 


5 Numerical simulations and discussions 

In the previous model [26], we have shown that the use of a vaccine with efficacy 
of about 60%, was to be accompanied by other measurements control such as 
means of personal protection (Information in relation to the damage caused by 
these diseases, spanning wearing clothes during hours of vector activity, use of 
repellents), vector control (combinig the use of Adulticide to kill adult vectors, 
chemical control with use of Larvicide to kill the eggs and larvae, and mechanical 
control to reduce the number of breeding sites at least near inhabited areas) [15]. 
Here, we investigate and compare numerical results, with the different scenario. 
We use the following initial state variables ^^(O) = 700, 14(0) = 10, Eh{0) = 220, 
4(0) = 100, RhiO) = 60, 4(0) = 3000, 4(0) = 400, 4(0) = 120, E(0) = 10000, 
L = 5000, P = 3000. 
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Table 6: Partial Rank Correlation Coefficients between Rq and each parameters of 
model (3). 


Parameter 

Correlation 

Coefficients 

Parameter 

Correlation 

Coefficients 

Parameter 

Correlation 

Coefficients 

tti 

-0.6125 

1 

0.3767 

Iv 

0.0378 

Oi2 

0.5960 

e 

-0.3348 

f^L 

-0.0357 

(^hv 

0.5817 

s 

0.2945 

Cm 

-0.0271 

Pvh 

0.5815 

a 

-0.1808 

Vh 

0.0178 

9 

0.5078 


-0.1594 

Vi 

-0.0161 

a 

0.4810 


0.1306 

fJ'E 

-0.0113 

fJjy 

-0.3911 

Ih 

-0.0605 

e 

-0.0109 

Tl 

0.4195 

Vv 

0.0578 

5 

-0.0077 

Te 

0.3888 


0.0439 

V2 

0.0037 

Afc 

-0.3876 

UJ 

0.0410 




5.1 Strategy A: Vaccination combined with individual pro¬ 
tection only 

In this strategy, we consider the model (3) withont vector control, we set 02 = 1 
and Cm = Vi = V 2 = 0 and vary the parameter related to individual potection, 
namely ai, between 0 and 0.8. The valnes of other parameters are given in Table 
4. Fignre 8 shows that the increase of the individnal protection level, permit to 
reduce the total number of infected humans, and the total number of infected 
vectors, but has no impact on the populations of eggs and larvae. However, from 
this hgure, it is clear that, this reduction is signihcant if the level of protection must 
turn around 80% at least, and this, over a long period. Thus, continuous education 
campaigns of people, on how to protect themselves individually, are important in 
the hght against the spread of arboviral diseases. 

5.2 Strategy B: Vaccination combined with adulticide 

Nowadays, Deltamethrin is the most used insecticide for impregnation of bednets, 
because it is a highly effective compound on mosquitoes at of very low doses [55]. 
However, when sprayed in an open environment, Deltamethrin seems to be effec¬ 
tive only during a couple of hours [56, 15]. Also, its use over a long period and 
continuously, leads to strong resistance of the wild populations of Aedes aegypti. 


36 





a 




gamma h 
eta_v 


gai-■■ 



-1 - 0.8 - 0.6 - 0.4 - 0.2 


0 


0.2 


0.4 


0.6 


0.8 


1 


(A) 


alpha_1 


alpha_2 

beta_hv 

beta_vh 

theta 

a 

mu_v 

Gamma_L 

Gamma_E 

Lambda_h 

I 

epsilon 

s 

sigma 

mu_P 

mu_h 

gamma_h 

eta_v 

mu_b 

omega 

gamma_v 

mu_L 

c_m 



eta_h 

eta_1 

mu_E 


xi 


delta 

eta_2 


- 0.8 


- 0.6 


- 0.4 


- 0.2 


0 


0.2 


0.4 


0.6 


0.8 


(B) 


Figure 7: Local (A) and global (B) sensitivity indices for Rq against model param¬ 
eters show that the local sensitivity results are not robust: the order of the most 
important parameters for Rq from the local sensitivity analysis not match those 
from the global sensitivity analysis. 



20% and 80%. To be more realistic, we will consider the technique called ’’pulse 
control” (the control is not continuous in time order is effective only one day every 
T days) [15]. To this aims, we consider that spraying is carried out once a week. 
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Figure 8: Simulations results showing how the total number of infected humans 
and the total number of infected vectors decrease when the individual protection 
increase. All others parameters values are in Table 4. 


and this, for 100 days. We set ai = rji = r ]2 = 0 and 02 = 1. 

Simulation result on hgure 9 show that a mortality rate induced by the use of 
larvicide, cm, greater than 60% has a signihcant impact on the decrease of the total 
number of infected humans and vectors, and on the decrease of eggs and larvae. 
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Figure 9: Simulations results showing how the total number of infected vectors, 
eggs and larvae populations dicrease when the aldulticide control parameter Cm 
increase. All others parameters values are in Table 4. 


5.3 Strategy C: Vaccination combined with larvicide 

Since the efficacy and the duration of a larvicide {Bti=Bacillus thuringiensis var. 
israelensis) strongly depend on several factors like water quality, exposure, and 
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even the type of breeding sites. To be more realistic, we thus consider that the 
duration can vary between a couple of days and two weeks [15, 57]. We consider 
that the larvicide spraying happens once every 15 days, and this, on a period of 
100 days. We set ai = Cm = 0 and 0^2 = 1- 

The hgure 10 shows that the use of larvicide has no signihcant impact on 
the decrease of total number of infected humans and vectors, as well as on the 
number of eggs and larvae. This can be justihed by the fact that the use of 
conventional larvicides neccssite certain constraints on their use: they can not be 
used continuously, their duration of action decreases with time. In addition, eggs 
of certain populations of vectors such as Aedes albopictus, come into prolonged 
hibernation when conditions in the breading sites are not conducive to their good 
growth (this is justified by the control rate value rji = 0.001). Also, the pupae do 
not consume anything, until reaching the mature stage. 

5.4 Strategy D: Vaccination combined with mechanical con¬ 
trol 

The effectiveness of this type of control depends largely of the awareness campaigns 
of local populations, in the sense that, to reduce the proliferation of vectors, people 
must always keep their environment clean by the systematic destruction of breeding 
sites. So, we consider that this type of control can be achieved by local populations, 
and this, every daily. We set ai = Cm = 0 = r/i = 772 - 

The hgure 11 shows that this type of control is appropriate in the hght against 
the proliferation of vectors. This can only be possible by the multiplication of local 
populations awareness campaigns. 

5.5 Strategy E: Combining vaccination, individual protec¬ 
tion and adulticide 

In this strategy, we consider the model (3) without larvicide and mechanical con¬ 
trol. we set a 2 = 1 and 7/1 = 7/2 = 0 and vary the parameter related to individual 
potection and the use of adulticide, namely ai and c^, respectively, between 0 and 
0.8. The values of other parameters are given in Table 4. Figure 12 shows that the 
use of the combination of these controls decreases signihcantly the total number of 
infected humans, infected vectors as well as the number of eggs and larvae, when 
its associated rates, namely ai and Cm, are greater than 0.3 and 0.2, respectively. 
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Figure 10: Simulations results showing how the total number of infected huamans, 
the total number of infected vectors, and the eggs and larvae populations dicrease 
whith the larvicide control associated parameters iji and 772 . All others parameters 
values are in Table 4. 
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Figure 11: Simulations results showing how the total number of infected vectors, 
eggs and larvae populations dicrease whith the mechanical control associated pa¬ 
rameter a 2 - All others parameters values are in Table 4. 


5.6 Strategy F: Combining vaccination, individual protec¬ 
tion and mechanical control 

Like for strategy E, the combined use of these three types of controls has a positive 
impact in the vector control. 
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Figure 12: Simulations results showing the advantage that we have to combine 
vaccination, individual protection and adulticide. 


6 Conclusion 

In this paper, we derived and analyzed a deterministic model for the transmission of 
arboviral diseases with non linear form of infection and complete stage structured 
model for vectors, which takes into account a vaccination with waning immunity. 
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Figure 13: Simulations results showing the advantage that we have to combine 
vaccination, individual protection and mechanical control. 


individual protection and vector control strategies. 

We begin by calculated the net reproductive number M and the basic repro¬ 
duction number, Rq, and investigated the existence and stability of equilibria. The 
stability analysis reveals that for A/" < 1, the trivial equilibrium is globally asymp¬ 
totically stable. When AT > 1 and i?o < 1, the disease-free equilibrium is locally 
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asymptotically stable. Under certain condition, the disease-free equilibrium is also 
globally asymptotically stable. We found that the model exhibits backward bifur¬ 
cation. The epidemiological implication of this phenomenon is that for effective 
eradication and control of diseases, Rq should be less than a critical values less 
than one. Thus, we proved, that the disease-induced death is the principal cause 
of the backward bifurcation phenomenon, in the full model and the corresponding 
model without vaccination. However, the substitution of standard incidence with 
mass action incidence removes the backward bifurcation phenomenon. 

We proved that the model admits at least one endemic equilibrium, and only 
one endemic equilibrium point in the model without disease-induced death, and in 
the model with mass action incidences, whenever the basic reproduction number 
is great than unity. 

Using parameters value of Chikungunya and Dengue fever, we calculated the 
sensitivity indices of the basic reproduction number, i?o, to the parameters in the 
model using both local and global methods. Local sensitivity analysis showed that 
the model system is most sensitive to a, the average number of mosquitoes bites, 
followed by the natural mortality rate of vectors. Considering that all input 
parameters vary simultaneously, we use the Latin Hypercube Sampling (LHS) to 
estimate statistically the mean value of the basic reproduction number. The result 
showed that the model is in an endemic state, since the mean of Rq is 2.0642, 
which is greater than unity. Then, using global sensitivity analyisis, we computed 
the Partial Rank Correlation Coefficients between Rq and each parameter of the 
model. Unlike the local sensitivity analysis, the global analysis showed that the 
parameters ai, the human protection rate, has the highest influence on Rq. The 
other parameter with an important effect are 0 : 2 , the efficacity of the mechanical 
control, Phv, the probability of transmission of infection from an infectious human 
to a susceptible vector, 13^^, the probability of transmission of infection from an 
infectious vector to a susceptible human, and 6 , the maturation rate from pupae 
to adult vectors. This showed that the order of the most important parameters for 
Ro from the local sensitivity analysis not match those from the global sensitivity 
analysis. So, the local sensitivity results are not robust. 

To assess the impact of combination of different controls, we conduct several 
simulations, using the called ’’pulse control” technique. According to the numerical 
results, we conclude that the use of an imperfect vaccine with low efficiency com¬ 
bined with high individual protection and good vector control statgie (reduction of 
breeding sites by local populations action, chemical action and use of adulticide), 
can effectively reduce the transmission of the pathogen and the proliferation of 
vector populations. However, due to lack of resources to implement these control 
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mechanisms, developing countries should focus on the education of the local popu¬ 
lations. Because, unlike diseases such as malaria whose breeding sites of Anopheles 
mosquitoes are known, those of arboviruses (old tires, flower pots, vases and other 
hollow...) are smaller and unknown for many local populations, which favor the 
development of vectors. 

Thus, pending the development of a high efficacy vaccine and long-acting, in¬ 
dividual protection and the various vector control methods are effective ways to 
overcome the arboviruses, for developing countries. In addition, the realization of 
these combination of controls may be too expensive, because it means that, for 
constant controls, we must keep them at levels high, and this, for a long time. 
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A Usefull result. 

We use the following result to compute the threshold at Eq. (11). 

Lemma A.l ([37]). Let M he a square Metzler matrix written in block form 
C D ) ^ ^ square matrices. M is Metzler stable if and only if 

matrices A and D — CA~^B are Metzler stable. 


B Proof of Theorem 3.1. 


The Jacobian matrix of / at the Trivial equilibrium is given by 



(27) 
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where 


Dh = 


( -^1 
e 
0 
0 
0 
0 


UJ 

-k2 

0 

0 

0 

0 


0 

0 

-h 

Ih 

0 

0 


0 

0 

0 

—A;4 

a 

0 


0 

0 

0 

0 


0 

0 

0 

0 

0 




Df2 = 


0 

~ / 

a(l - ai)/3hvS^ 
a(l - ai)'$hv'^ri^V^ a(l - ai)/3ft„7rV'^° 


V 

/ a(l - ai)f3hvVvSh 



( 0 

0 

0 

0 

0 

0 

\ 


0 

0 

0 

0 

0 

0 


,Df3 = 
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0 

0 

0 

0 

/ifc 
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0 

0 

0 

0 

0 



VO 

0 

0 

0 

0 

0 

/ 


0 0 0 








a(l - ai)$hvVvH° o(l - ai)/3hvH^ 


0 

0 

0 


0 

0 

0 


0 0 0 


0 0 0 


0 0 
0 0 


0 


Dh = 


/ 


V 


-kg 

Iv 

hfe 

0 

0 


0 

-k^ 

fJ'b 

0 

0 


0 

0 

-^5 

s 

0 


0 

0 

0 

-ke 

I 


0 \ 
0 
0 
0 

-kr / 


0 
0 

0 9/ 


and = ^0 + ttK?. 


The characteristic polynomial of Df(So) is given by: 

P(A) = — (A + ks) (A + ^ 4 ) (A + kg) (A + kg) (A + Hh) 0i(A)02(A) 


where 

0i(A) = A^ + {k 2 + ki)X + Hh{k 2 + 0 and 02(A) 

we have set 



^1 — k^ -\- kg -\- ki -\- kg, A2 — k^{k^ + fcg + k^) + k^{k^ + kg) + k^kg, 

A 3 = k^kgkj + k^ik^k^ + k^ik^ + /cg)), A^ = k 3 k^k'jk 3 {l — J\f). 


The roots of P(A) are Ai = —/ih, Ai = —ki, X 2 = —ks, A 3 = —k^, X 4 = —kg, 
A 4 = —kg, and the others roots are the roots of 0i(A) and 02 (A). The real part 
of roots of 01 (A) are negative. Since A/” < 1, it is clear that all coefficients of 
02 (A) are always positive. Now we just have to verify that the Routh-Hurwitz 


criterion holds for polynomial 02 (A). To this aim, setting Hi = Ai, H 2 


Ai 1 
A3 A2 ’ 
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1 0 0 

As A 2 Ai 1 

0 A 4 A 2 

0 0 0 A4 

The Routh-Hurwitz criterion of stability of the trivial equilibrium is given by 



Ai 

1 

0 


Hg = 

As 

A2 

Ai 

,H 4 = 


0 

A4 

CO 



{ ifi > 0 f Hi >0 

H 2 > 0 I H 2 > 0 

if 3 > 0 iis > 0 

if4 > 0 [ y44 > 0 


(28) 


We have Hi = Ai = + kg > 0, 


H 2 = A 1 A 2 — As 

= (/c7 + /cg + fcg) kg + (fcy + (2 /c 6 + 2/C5) /c7 + fcg + 2 k^kQ + /cg) /cg 
+ (fcg + fcs) ^7 + (fcg + 2k^kQ + fcg) A;7 + k^k^ + k‘^k% 


H 3 = TiR2^3 - AIA 4 - Al 


= (fee + h) {k^ + (/C6 + kA A + k^ko) A:| 

+ (^fXb^sO + (fcg + kA kj + 2(/c6 + kA^^k^ + (fcg + 4^5^^ + 4 /C 5 /C 6 4- /c^) A ~k k^k^ 4- ^k^k^ k^kA) kg 

4" [(2/^7 4“ 2kQ 4“ 2kA fiAsO 4“ (fcg 4" ^k^kg -\- fc^) k^ 4“ (A^g 4“ 4/^5^^ 4“ Ak‘3kQ -\- fcg) 

4“ (Ak^kg 4“ Akgkg -\- 2k^kA) A 4“ kgk^ -\- /cg/cg] fcg 4“ {k'j -\- (2/cg + 2kA k'^ -\- kg A ^k^kQ -\- kA) fib^s9 
+ {kbkl 4- fcg/cg) ky 4- (fcgfcl 4- 2 / 05 ^! 4- fcf/cg) k^ + {k^kl 4- /cf/cg] ^7 

We always have Hi > 0, H 2 > 0, Hg > 0 and ii 4 > 0 if Ai < 1. Thus, the trivial 

equilibrium £q is locally asymptotically stable whenever Ai < 1. 

We assume the net reproductive number Ai > 1. Following the procedure and 
the notation in [35], we may obtain the basic reproduction number Rq as the dom¬ 
inant eigenvalue of the next-generation matrix [34, 35]. Observe that model (3) 
has four infected populations, namely E^, R, E^, R- If follows that the matrices 
F and V defined in [35], which take into account the new infection terms and 
remaining transfer terms, respectively, are given by 

/ a(l - aAAhvVvH^ a{l - aAAhvH^ \ 


000 

a(l - aAAvhVvS^ a(l - aARuS^ 

K K 

\ 0 0 


0 

0 


/ 
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0 













( h 0 0 0 \ 

^ - 7 fe ^4 0 0 

0 0 fcg 0 ■ 

\ 0 0 - 7 ^ ks 

The dominant eigenvalue of the next-generation matrix FV~^ is given by (8). 
The local stability of the disease-free equilibrium £i is a direct consequence of 
Theorem 2 of [35]. This ends the proof. 


C Proof of Theorem 3.2. 


Setting Y = X -TE with X = {Sh,Vh, E^, h, Rk, S,, E„ I,, E, L, Pf, = 


(7 + nV«), 4* = 
rewrite (3) in the following manner 


E 

, and Aio = ( ^6 + s— 

El 


we can 


dY 

dt 


B{Y)Y 


(29) 


where B{Y) 


A{Y) B{Y) 
C{Y) D{Y) 
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It is clear that V = (0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) is the only equilibrium. 

Then it suffices to consider the following Lyapunov function £(V) =< g,Y > 


were g = 


1 , 1 , 1 , 1 , 1 , 1 , 1 , 1 , 


fib flbSl 


Straightforward computations 


lead that 


t{Y) =<g,Y>‘'^<g,B{Y)Y> 

P-Zihi fihY2 g^hY^ (h/i T YfY^ giiYb 


-^{Y, + Yj + Ys) 

Ke 


hks 

gbKi 


Igddo + 0 



We have C{Y) < 0 if < 1 and C{Y) = 0 if F* = 0, i = 1, 2,..., 11 (i.e Sh = 

Vh = and Eh = Ih = Rh = Sy = Ey = ly = E = L = P = 0). Moreover, 

the maximal invariant set contained in |£|£(y) = o| is (0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0). 

Thus, from Lyapunov theory, we deduce that (0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0) and thus. 
So, is GAS if and only J\f < 1. 


D Proof of Theorem 3.5. 


In order to determine the existence of endemic equilibria, i. e. equilibria with all 
positive components, say 


_ / o* T/'* IP* r* D* o* IP* r* rp r d\ 

^ ^ ^ 


we have to look for the solution of the algebraic system of equations obtained by 
equating the right sides of system (3) to zero. In this way we consider two case: 


(i) Special case: Absence of disease-induced death in human (5 = 0 ) 

Note that in the absence of disease-induced death in human population, we have 
^h = ^h= Rh/gh- Let 

^C,* _ a(l - ai)/3hy{VyE; + /*) _ a(l - ai)l3yh{r]hEl + ID 

K - -“- Ni - 
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be the forces of infection of humans and vectors at steady state, respectively. 
Solving the equations in (3) at steady state gives 


O. A,(vrAr + k,) jSl 

^^h{k2 + + XjI +7r/ci + fc2) ^ +^ 2 ) 


El = 


>^r{Sl + TiV^) 


T* — 


hh 


ID* _ 

— 


a-fhXh*iE*h + T^Vh. 


and 


9P 

{Xr + ksY 


E = 


^bOKEP 


_epx^ 
" kYxr + ksY 
L = 


_ jjpxr 

kskg{Xv* + ksY 
libOsKEKiP 


(31) 


(32) 


{k^ksKE + iibOP) ’ keKiikbksKE + fib9P) + subOKEP' 
where P is solution of the following equation 

/(F) = -krP [i^b9{sKE + k^KL)P + k^k^ksKEKLiAf - 1)] = 0 (33) 

A direct resolution of the above equation give F = 0 or F = . 

IJ,b9{sKE + kaKL) 

Note that F = 0 corresponds to the trivial equilibrium Sq. Now we consider 
F > 0 i.e. M > 1. Replacing (31) and (32) in (30) give 


\C,* _ 

\ - 


( 2(1 OiYYhvf^h 

Ah, 


Vt 


9px: 


+ 


7v9px: 


kgiXl + kg) kskg{Xl + ks) 


xr = 


a{l-ai)Yvhf^h XliSl + TiV,:) , 7hA^(F^ + 7rR, 


A/] 


Vh 


k?. 


+ 




3 "^4 


Substuting (35) in (34) give 

{k^Ki + sKe) Ah [a2(Ah)^ + ^lA^ + Oo] = 0 
where 02, ai and oq are given by 

(tT^ + k2) ( (fclTT + /C2) a(l - OlYYvhYlh + kiT]h){TT^ + ^2) 


(34) 

(35) 

(36) 


Rh — 


7r(^ + ^ 2 ) 


f^h 


+ 


ksk^ks 


02 = (a(l - aYYvhl^hhh + kir]h) + k^k^kYj fcgpfeAhTT, 
k^k^kgkg^bX-h (^ F kY) r ^ d ^ 

- \J^b — -Kij, 


( 37 ) 


(1\ — 


(tT^ + kg) 

ao = iihk^kAkskgiibXh{i + kg) (1 - Fi) 
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The trivial solution = 0 of (36) corresponds to the disease-free equilibrium 
E\. Now, we just look the equilibria when > 0. Note that coefficient 02 is 
always positive and uq is less (resp. greather) than unity if and only if i?i > 1 
(resp. R\ < 1). Thus model system (3), in absence of disease-induced death in 
human population (5 = 0), admits only one endemic equilibrium whenever Rq > 1. 
Since the sign of coefficient Oi depend of the value of parameter, we investigate 
the possibility of occurence of backward bifurcation phenomenon when i?o < 1 - 
Furthermore, consider the inequality 

Ri El Rb- ( 38 ) 


Since 02 is always positive and oq is always positive whenever Rq < then, the 
occurence of backward bifurcation phenomenon depend of the sign of coefficient oi. 
The coefficent oi is always positive if and only if condition (38) holds (i.e Ri < Rh). 
It follows that the disease-free equilibrium is the unique equilibrium when J\f > 1 
and Rq < 1. Now if R^ < Ri < 1, then in addition to the DFE 81 , there exists 
two endemic equilibria whenever A = af — 4 a 2 ao > 0. However, Rb < Ri < 


Rb < 1 ^ < 


+ Hh) tt) ^ + (cu + fihf] ksk^kg 


< 0 . 


a(l - ai)/ife(7r^ + ^ 2 ) 2 ( 7 /* + k^iqh) 

Since all parameter of model (3) are nonnegative, we conclude that the condition 
Rb < Ri < 1 does not hold. And thus, the backward bifurcation never occurs in 
the absence of disease-induced death in human. 


(ii) Presence of disease induced death in human (5 7 ^ 0). In this case, we have 
. Applying the same procedure as case (i), we obtain that A^ at 


Nl = 


- ^Ih 

l^h 


steady state is solution of the following equation 

/(-^D = K [c4(A)))'^ + C 3 (A)))^ + C 2 (A’^)^ + CiA^ + Co] = 0, 


(39) 


where 


C 4 = kgK 12 fib^h {kzki - Sjh) {kioafXhil - ai)f3yh + ks{ksk4^ - 5'jh )), 

C 3 = TT^kskik^kekiokiia^- aif jShvmrPvhKEKL + 2kQkioKi2afj.bdAhfJ.hlhTr{^ - ai)j3yhi 

- kzkikgkiQKi2a^bAhiJihT^{^ - ai)l3vh^ - ‘2kskgKi2fj,bS‘^Ah"flTT^ + 2k3k4kskgKi2HbSAhjhTr^ 

- kik3k4kgkioKi2aHbAhfJ.h'^{^ - ai)/3«/i + ‘^k2kgkioKi2aHbSAbHhlh(,^ - ai)l3vh 

- 2k2k3k4k9kioKi2af4bAhfJ.h{^ - a:i)l3yh + 2kik3k4k3kgKi2^ih&Ah'yhT^ - 2kiklklk3kgKi2iibAhTT 

- 2k2kskgKi2HbS‘^Ahjl + 4k2k3k4kskgKi2fib5Ahjh - 2k2klklk3kgKl2^ibAh), 
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C2 = kiki^k^kek^okiia^- ai)^l^vh^K eK 

+ kikskik^k^kiokiia^ fij^{l — ai)^ jS^f^K eK e 

+ 2k2kzkikzkQkiQkiia'^ ^i\{l — ai)'^ eK l + kQkiQKi20-IJ-b^^hk-hlh'’^'^ — ch)/ 3 d)i?^ 

— kskQKi2tJ.bS^Ah'ylTT^^^ - kik3k4kQkioKi2atif,A,^tii^E^{l - 

+ k3k4kgkioKi2atJ.b^hlJ'h‘^'^{^ - ai)l^vh^ + 2k2kgkioK42at4biAhlJ.h'yhTr{i - ai)l3^f,i 

— k2k3k4kgkigK42aiJ.bf^hk-h'^(^ - 0‘l)Pvhi + 2kik3k4k3kgKi2t4biAh-yh'^^^ 

— 2k3k4k3kgKi2^i.b6^hlhb^'’^i + 2k3k\k3kgK42f^b^h'^'’^^ — '^kgkskgKigfi-b^^ ^hAh'^'S 

+ 2k2k3k4kskgK42lJ'bSAh,'yhTr^ - 2kik2k3k4kgk4(}Ki2aiibAf,t^h{^ - ctl)'^Pvh 

+ kgkgkioKigaiibSAht^h'yhi^ — cil)0vh - k2k3k4kgkioKi2afibAb,tih{^ — cil)0vh 

— kfk'^k^kskgKigf^b^h'’^'^ + 4A:ifc2fe3fc4fc8fc9-f4’l2Ai6<5A^7/i7r — ikikgk'^k^kskgKig^ib^h'’^ 

— kgk^kgKigfJ'bS^ + ‘^f^2^3k4kskgKi2lJ.bSAb,'yh — fc2^3^4^8fc9Ari2At(,A/i, 

Cl = {(kik3k4k5kekiokiia‘^ - ai)l3hvmr‘^ + ksk^k^k^kiokiia? ijl\{ 1 - aifl3hvn{k2 - uj)Tr)l3yhi 

+ {2kik2k3k4k5kekiokiia‘^- ai)l3hvmv + klk3k4k3k3kiokiia^y,\{l - ai)/3?i„n)(l - ai)l3yh)KEKL 

+ {k3k4kgkioKi2aHb^hk-h^T^{^ - ai)Pvh - 2k3k4kskgKi2fJ,bSAh'yh^'^)^^ 

+ {{k2k3k4kgkiQKi2afibAhtJ-hi^ - kik2k3k4kgkioKi2afj,bAhHhE)il - ai)/3vh 
+ {2kiklklkskgKi2fJ,bAhUj + 2kik2k3k4kskgKi2fibSAhjh)E 
+ { 2 k 2 klklkskgKi 2 fibAh - 2k2k3k4kskgKi2fJ.bSAh^h)^^)^ 

- kiklk3k4kgkioKi2a^ibAh^ih(.^ - ai)l3vh - 2klk2klklk8kgKi2fibAhTT 
+ 2kik^k3k4kskgKi2fibSAh"fh - 2 kik^klklkskgKi 2 fJ.bAh, 


Co = klklkskgKi2iJ,b^hlJ^l{k2 + 0^ (^o “ l) ) 

with fcio = 7h + Vhki, fell = 7i, + Vvks, = {sKe + kgRi) and n = M -1. Notes 
that C 4 is always negative and Cq is positive (resp. negative) if Rg is greather (resp. 
less) that the nnity. It follows, depending of the sign of coefficients C 3 , C 2 and 
Cl, that the model system (3) admits at least one endemic eqnilibrinm whenever 
i?o > 1 and the phenomenon of backward (resp. forward) bifurcation can occurs 
when i?o < 1 (resp. Rq > 1). This ends the proof. 
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